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abstract 

We  present  a  population  of  56  linear,  two-dimensional  elliptic  partial 
differential  equations  (PDEs)  suitable  for  evaluating  numerical  methods  and 
software.  Forty  two  of  the  PDEs  are  parameterized  which  allows  much  larger 
populations  to  bemade;  189  specific  cases  are  presented  here  along  with  solu¬ 
tions  (some  are  only  approximate) .  Many  of  the  PDEs  are  artifically  created 
so  as  to  exhibit  various  mathematical  behaviors  of  interest;  the  others  are 
taken  from  "real  world"  problems  in  various  ways.  The  population  has  been 
structured  by  introducing  measures  of  complexity  of  the  operator,  boundary 
conditions,  solution  and  problem.  The  PDEs  are  first  presented  in  mathematical 
terms  along  with  contour  plots  of  the  189  specific  solutions.  Machine  readable 
descriptions  are  given  in  Part  2,  I  IRC  Technical  Summary  Report  #2079;  many  of 
the  PDEs  involve  lengthy  expressions  and  about  a  dozen  involve  extensive 
tabulations  of  approximate  solutions. 
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1.  INTRODUCTION 


The  motivation  for  creating  this  PDE  population  is  for  use  in  the  evaluation 


of  numerical  methods  and  PDE  software.  The  need  and  rationale  for  a  systematic 


approach  to  such  evaluations  is  given  in  (Rice,  1979],  (Houstis  and  Rice,  1980], 


[Crowder,  Dembo  and  Mulvey,  1979];  it  suffices  here  to  say  that  a  properly  chosen 


problem  population  is  an  essential  ingredient  for  a  sound  evaluation  of  numerical 


methods  and  software. 


A  useful  population  of  PDEs  is  inevitably  very  lengthy  and  this  one  is  no 


exception  as  one  sees  from  the  last  two  appendices.  Thus  in  the  body  of  this  paper 


we  discuss  the  sources  of  the  PDEs ,  how  they  are  described  in  the  appendices  and 


how  a  structure  has  been  created  in  the  population  through  the  use  of  quantitative 


(but  subjective)  measures  of  features. 


It  is  important  that  one  be  able  to  create  relevant  subpopulations  as  one 


inevitably  wants  to  evaluate  methods  for  particular  subclasses  of  PDEs  (e.g.. 


separable,  with  singularities  or  with  mixed  boundary  conditions).  Experience  shows 


that  no  one  universal  method  is  nest  for  all  PDEs  (even  in  this  rather  restricted 


context)  and  one  of  the  important  tasks  of  research  is  to  create  and/or  identify 


methods  that  are  especially  efficient  for  particular  classes  of  PDEs.  Once  one 


embarks  on  such  a  task  one  sees  that  this  population,  which  originally  might  seem 


large  and  bulky,  is  actually  rather  small  for  the  uses  to  be  made  of  it.  It  is 


only  the  fact  that  it  can  be  substantially  expanded  in  various  directions  through 


the  parameterization  that  gives  one  hope  that  it  is  adequate  for  a  wide  variety  of 


evaluations . 


Sponsored  by  ti.e  United  States  Army  under  Contract  Nos.  DAAG29-75-C-0024  and 
DAAG29-8Q-C-0J41.  This  material  is  based  upon  work  supported  by  the  National 
Science  Foundation  under  Grant  Nos.  MCS77-01403  and  MCS79-01437. 
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2.  characteristics  of  the  problems. 

A  source  parameter  is  assigned  to  each  PDE  which  ranges  from  0  (artificial 


problem)  to  100  (actual  real  world  problem) .  This  feature,  as  the  others  introduced 


later,  is  subjective  in  nature  and  the  values  given  must  be  taken  as  approximate 

indications  of  our  intuitive  feelings.  The  PDE  u  +u  =1  might  be  completely 
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artificial  for  one  person  and  be  the  actual  applications  PDE  for  another.  We  have 
at  least  tried  to  be  consistant  in  these  values. 

2.1.  Sources 

Many  problems  have  been  normalized  so  the  maximum  value  of  the  solution  is 
1.0  and  all  have  this  value  between  .1  and  100.  Many  of  the  domains  have  been 
standardized  to  the  unit  square,  0  <_  x,y  1  .  The  sources  of  the  PDEs  are: 

A.  Problems  used  in  previous  studies.  Nine  problems  are  included  which 
were  used  by  [Eisenstat  and  Schultz,  1973]  or  [Houstis  et  al,  1975  and  1978]. 
Subsets  of  this  population  have  been  used  by  [Houstis  and  Papatheodorou,  1977  and 
1979]  and  [Lynch  and  Rice,  1978],  Some  of  these  PDEs  have  had  parameters  added 
and  all  have  been  normalized  so  the  maximum  value  of  the  solution  is  about  1.0. 

B.  Artifical  Problems.  Many  problems  have  been  created  just  to  exhibit 
various  mathematical  behaviors  of  interest  (e.g.  singularities,  oscillations  or 
wave  fronts).  Such  behaviors  are  important  for  theory  or  application  (or  both) 
and  one  needs  to  have  them  present  in  the  population  in  an  easily  identifiable 
manner . 

C.  Problems  adapted  from  the  "real  world",  a  persistent  difficulty  is  the 
desire  to  have  PDEs  which  represent  the  "real  world"  and  the  necessity  to  know 
their  true  solutions.  Among  the  strategies  to  adapt  real  'world  problems  we 

have  used: 

(i)  choosing  explicit  functions  which  model  the  physical  solutions  and  the: 


determining  appropriate  boundary  conditions  and/or  right  side  to  make 
this  the  true  solution. 


{ ii )  using  truncated  series  expansions  (of  high  accuracy)  with  appropriate 
small  modifications  in  the  boundary  conditions  or  right  side. 

(iii)  solving  nonlinear  problems  approximately,  then  substituting  the 

tabulated  numerical  solution  into  the  operator  (using  quadratic  inter¬ 
polation  from  a  10  by  10  grid)  to  obtain  a  linear  problem  which  is,  in 
turn,  solved  approximately.  In  these  cases  the  true  solution  is  not 
known,  but  the  machine  readable  population  contains  tabulated  values  of 
a  hopefully  accurate  numerical  solution. 

2.2  Problem  Features  and  Complexity  Classifications.  We  identify  as  problem 
features  the  smoothness  and  local  variation  of  operator,  the  boundary’  conditions  and 
the  solution.  These  features  are  quantified  on  a  one-dimensional  scale  of  0  to  I  TO 
even  though  there  are  rather  independent  properties  that  can  be  called  smoothness  or 
local  variation.  These  features  are  measured  subjectively  from  the  following  de¬ 
scriptions  of  the  scale. 

Smoothness .  This  refers  to  the  mathematical  properties  of  the  functions  or 
operators  involved.  Key  points  on  the  scale  are: 

00  =  entire  functions  or  constants 
1C  =  analytic;  very  well  behaved 

30  =  very  smooth,  some  higher  derivative  (5  or  so)  discontinuity  possible 
50  =  still  smooth,  third  derivative  discontinuity/  possible 
70  =  not  rough  to  the  eye,  but  possibly  only’  1  continuous  derivative 
80  =  continuous,  functions  might  be  theoretically  smooth  but  rough  on 
a  gross  scale 

90  =  possibly  discontinuous,  nearly  singular  functions  or  operators 

2 

100  =  strong  singularities  like  1/x  or  1/x"  . 

Local  variation.  This  refers  to  how  much  a  function  changes  (relative  to 
its  size)  in  a  small  part  of  its  domain.  These  variations  might  be  oscillations, 
wave  fronts,  peaks  or  boundary  layers.  Key’  points  on  the  scale  are: 
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00  =  very  smooth,  uniform 

10  =  mild  variation,  probably  convex,  some  non-uniformity,  e.g. 
sin(2x),  e3x  on  [0,1] 

25  =  modest  variation  of  oscillation;  mild  wave  front  or  peak,  e.g. 
sin(6x),  l/(l+100x^)  on  [0,1] 

40  =  considerable  peak  or  oscillation;  change  of  magnitude  occurs  within 
10-15%  of  domain 

60  =  sharp  peaks,  wave  fronts,  boundary  layers  or  oscillations;  100%  change 
in  magnitude  occurs  within  5%  of  domain 

75  =  practically  a  discontinuity  in  magnitude;  continuity  observable  only 
with  a  fine  scale  examination 

90  =  actual  discontinuity  in  magnitide;  extreme  oscillation,  step  functions, 
e.g.  SIN (300x)  on  [0,1] 

The  overall  problem  complexity  is  represented  by  the  average  of  the  above 
six  feature  measures.  The  problems  in  this  population  do  not  have  complexities 
exceeding  58  (only  one  exceeds  50) ,  a  level  which  might  be  interpreted  as  "rather 
messy  with  one  or  two  substantial  complications".  The  problem  feature  measures 
are  included  in  the  descriptions  along  with  the  source  parameter. 

Appendix  1  presents  some  summary  information  about  the  population.  Tables 
are  given  which 

A.  group  the  PDEs  according  to  types  of  the  operator  and  boundary  conditions 
(e.g.  Helmholtz  and  Dirichlet  or  constant  coefficients  and  mixed) 

B.  list  the  56  PDEs  with  abbreviated  feature  descriptions 

C.  group  the  PDEs  according  to  the  smoothness  of  the  operator  and  right  side 

D.  group  the  PDEs  according  to  the  smoothness  of  the  solution. 

A  figure  is  also  given  which  displays  the  overall  problem  complexity  for  the  189 


specific 
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3.  FORMAT  OF  PROBLEM  DESCRIPTIONS 

Appendix  2  contains  a  mathematical  description  of  each  PDE  along  with  contour 
plots  for  each  specific  instance  included  in  the  set  of  189  PDEs.  An  example  is 
shown  in  Figure  1.  The  description  begins  with  a  problem  number  and  source  followed 
by  a  matnematical  description  of  the  PDE.  Then  brief  comments  are  given  for  the 
operator,  right  side,  boundary  conditions,  solution  and  parameters  (if  any).  Some¬ 
times  functions  appearing  in  the  mathematical  description  are  defined  in  these 
comments . 

Four  generic  functions  are  used: 

f(x,y)  =  right  side  of  PDE  determined  so  that  the  given  true  solution  is 
correct. 

f (x) ,g (y)  =  right  sides  of  boundary  conditions  determined  so  that  the  given 
true  solution  is  correct. 

T(x,y)  =  the  true  solution,  used  in  the  coefficients  of  some  PDEs  derived 
from  nonlinear  problems. 

r(x,y)  =  an  approximate  solution  used  in  some  PDEs  whose  true  solution  is 
unknown . 

Contour  plots  are  given  for  one  or  more  particular  PDEs  for  each  problem. 

The  border  of  the  plots  contains  the  following  information: 

(i)  values  of  the  parameters 

(ii)  maximum  and  minimum  values  of  the  solution;  the  contours  are  equi- 
spaced  between  these  values. 

(iiil  the  classification  parameters  in  the  form 
S.P  01.02  B1.B2  S1.S2 

where 


p 

0l2 

:ator. 


a  *=  1.  10. 

ty  in  operator 
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Figure  1.  An  example  of  the  mathematical  description  of  a  PDE  along  with  some 

contours . 

The  machine  readable  description  of  the  PDE  population  consists  of  two 
files:  EQNFIL  and  MACFIL.  EQNFIL  has  189  entries  which  are  either  complete 
statements  of  the  ?DE  in  the  ELLPACK  language  (see  [Boisvert,  Houstis  and  Pice, 
1979])  or  a  reference  to  an  entry  in  MACFIL  with  values  given  for  parameters.  Se 
Figure  2  for  a  short  example.  The  information  -,iven  starts  with  the  problem  numb 
feature  parameter  values  and  a  code  for  various  attributes  of  the  PDE  which  are 
used  within  the  ELLPACK  system.  Then  ELLPACK  language  code  is  given  for  the  oper 
ator  and  boundary  conditions;  this  code  should  be  self  explanatory  once  one  sees 
that  UXX$  represents  u  ,  etc.  Finally,  there  is  a  Fortran  ''ode  for  any  function 
that  appear  in  the  operator,  right  side  or  boundary  conditions.  Thir-  latter  code 
averages  about  20  lines  and  can  be  as  much  as  150  lines  (excluding  tables  that  ar 
part  of  some  problems).  These  descriptions  are  giver,  in  Part  2  of  this  report. 
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MACFIL  entries  are  just  like  EQNFIL  descriptions  of  a  PDE  except  that  the 
places  where  parameter  values  are  to  be  substituted  are  indicated  by  SB,  etc. 

A  refers  to  the  first  parameter,  B  the  second  and  so  on.  There  are  somewhat  mor*. 
than  8500  lines  in  these  two  files. 


•  PROBLEM  2  • 


000.04  000.00  004. OS  010.02 

2000200000020 

TUO  DIMENSIONS 

UXXS  +  C1.*Y«YJUYYS  -  UX*  -  tl.+Y*YJUYS  *  FCX.YJ 
MIXED 

X=0.  .  MIXED  =  (l.)U  ♦  C  l.JUX  *  0.27*£XPCYJ 

X=l.  ,  MIXED  *  Cl.JU  ♦  C-l.JUX  =  0. 

Y=0.  .  MIXED  =  Cl.JU  ♦  C  l.JUY  =  0.27»EXPCXJ 

Y=l.  .  MIXED  =  Cl.JU  ♦  C-l.JUY  *  0.13S*CAL0GX2. J-l. J«CX»X-XJ«*2 

FUrCTIOM  TRUECX.YJ 

TRUE  *  0.13S*CEXPCX+YJ+CX«X-XJ**2«AL0GCl.'*Y»Y>i 

FETURM 

END 

FUNCTION  FCX.YJ 

F  *  0.135*C  C-4. •X»X»X+18. «X*X-14 .«X*2. J •flLOCC 1 -+Y«Y J 
*  -  2.«CCX»X-XJ*»2)»CY»Y+Y««3+Y-1.J^C1.+Y«Y)  J 

RETURN 
END 


•  PROBLEM  3  • 

•EOR 

•PARAHETER  SET  ICF^I.5) 

•  000.43  090. SO  000.00  070.40 

EXPflMD  3'1-S' 

•EOR 

•PARAttTTER  SET  2Cft=2.5J 

•  000.35  030.50  000.00  0S0.20 

EXPAND  tt.Sr 

•EOR 

•PARAMETER  SET  3Cft=3.5J 

•  000.28  070.30  000.00  050.15 

EXPAND  3^3.5/ 

•EOR 

•PARAMETER  SET  4(A=4.5J 

•  000.23  055.20  000.00  C40.20 

EXPAND  3^4.5^ _ _ 

Figure  2.  A  sample  from  EQNFIL  showing  a  short  PDE  description  in  machine  readt. ! 

form  and  a  reference  to  a  similar  description  in  MACFIL 
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APPENDIX  ONE:  TABULATIONS  OF  POPULATION  CHARACTERISTICS 


Table  1 

Classifications  of  Problems 

Accordinf  to  Operator  and  Boundary  Conditions 

Operator 

Constant  Coefficients 

Non-Constant  Coefficients 

Dirichlet 

Neumann 

Mixed 

Dirichlet 

Neumann 

Mixed 

Laplace 

3,  4,  7, 

B,  10,  11, 
17,  33,  34, 
35,  47,  50 

4,  31,  35, 
38.  55 

Helmholtz 

Type 

9.  41,  53 

6,  20,  39, 
44,  45,  48, 
49 

Self-Adjoint 

5 

1,  13,  22. 
25,  28,  54 

1,  19,  23. 
52 

General 

14,  46 

42 

43 

12,  15,  16, 
18,  21,  26, 
27,  29,  30, 
32,  36,  37, 
56 

_  -> 

*  - 

L - 

2,  23,  24, 
40,  51 

Note  that  problems  1,  4,  and  35  appear  in  two  places  in  the  table  since  they 
have  boundary  conditions  of  the  form 

u  +  auN  =  g 

and  hence  have  Dirichlet  boundary  conditions  for  a=0. 


Table  2 

Problem  Characteristics 


The  principal  characteristics  are  tabulated  below  using  the  following 
encodings: 


A  Analytic 
BL  Boundary  Layer 
C  Constant  (coefficients) 

CC  Computationally  Complex 
D  Dirichlet  Boundary  Condition 
E  Entire 
H  Homogeneous 
J  Jump  Discontinuity 
M  Mixed  Boundary  Condition 


N  Neumann  Boundary  Condition 
NS  Nearly  Singular 
0  Oscillatory 

P  Parameterized  or  Peaked 
S  Singular  (infinite) 

SD  Singular  Derivative 
U  Unknown 

VS  Variable  Smoothness 
WF  Wave  Front 


Problem 

Number 

Operator 

Right 

Side 

. 

Solution 

Boundary 

Conditions 

Domain 

IP 

A 

E 

E 

M 

Unit  Square 

2 

E 

A 

A 

M.H 

Unit  Square 

3P 

C 

S,SD 

S,SD 

D.H 

Unit  Square 

4P 

C 

E 

E 

M 

Unit  Square 

5P 

C 

E 

E 

D.H 

Unit  Square 

6 

E,NS 

A 

A.0 

D.H 

Unit  Square 

7 

C 

C 

SD 

D.H 

Unit  Square 

BP 

C 

SD 

SD.WF 

D 

Unit  Square 

9P 

C,NS 

E.NS 

E,  BL 

D 

Unit  Square 

10P 

C 

E.P 

E.P 

D.H 

Unit  Square 

IIP 

C 

A.0 

A.0 

D 

Unit  Square 

12P 

E,0 

E.O 

E.O 

D 

Unit  Square 

13 

J 

S 

SD 

D 

Unit  Square 

14P 

C 

S 

S 

D 

Unit  Square 

15P 

A,  NS 

S 

SD 

D 

Unit  Square 

18P 

A,  NS 

C 

U.BL 

D.H 

Variable  Square 

17P 

C 

A,  NS 

A.NS.WF 

D 

Unit  Square 

IBP 

E 

A,  NS 

A.NS.WF 

D 

Unit  Square 

19P 

S 

s 

E 

M.H 

Square 

20P 

NS.P.CC 

P 

E.P 

D 

Rectangle 

21 

E 

E 

E 

D 

Unit  Square 

22 

SD 

S 

E 

D 

Unit  Square 

23P 

SD 

SD 

SD.WF 

M.H 

Unit  Square 

24P 

S.NS 

S.NS 

U.P 

M.H 

Square 

25P 

SD 

S 

E 

D.H 

Unit  Square 

Table  2 

Problem  Characteristics 


Problem 

Number 

Operator 

Right 

Side 

Solution 

Boundary 

Conditions 

Domain 

26P 

A 

A 

U,  SD 

D.H 

Variable  Square 

27 

A.  NS 

C 

U.BL 

D.H 

Square 

28P 

J 

C 

U.WF 

D.H 

Square 

29P 

S 

H 

U.VS.BL 

D 

Unit  Square 

30P 

A.CC 

A,  CC 

A.  NS 

D 

Unit  Square 

31 

C 

C 

E.(SD) 

M 

Square 

32 

A 

A 

E 

D.H 

Rectangle 

33 

C 

E 

E.O 

D 

Rectangle 

34 

C 

C 

E.(SD) 

D 

Square 

35P 

C 

H 

E.O.BL 

V 

Square 

36P 

S 

S 

A.BL 

D 

Unit  Square 

37 

E 

E 

E 

D 

Unit  Square 

38P 

C 

H 

E.O.  VS 

D 

Rectangle 

39P 

CC.S 

CC.S 

U.BL 

D.C 

Unit  Square 

40P 

E 

A 

A 

M 

Unit  Square 

41P 

C,  NS 

SD.NS 

SD 

D.H 

Square 

42P 

C 

H 

A.O 

N 

Variable  Rectangle 

43 

C 

H 

E 

M 

Square 

44P 

CC 

CC 

U.BL 

D.H 

Unit  Square 

45P 

C,  NS 

H 

U.BL 

D 

Unit  Square 

46P 

C,NS 

H 

U.BL 

D 

Variable  Rectangle 

47P 

C 

S 

SD.VS 

D 

Unit  Square 

4BP 

CC 

CC 

U 

U 

Unit  Square 

49P 

CC 

CC 

U.SD.BL 

D.C 

Unit  Square 

50 

C 

H 

E.O 

D 

Rectangle 

51P 

S 

C 

U.SD.WF 

M.H 

Unit  Square 

52P 

CC 

H 

U.O 

M.C 

Unit  Square 

53P 

C.NS 

E.O 

E.O 

D 

Unit  Square 

54P 

E,  CC 

S.CC 

SD.VS 

D 

Unit  Square 

55P 

C 

H 

S.VS.BL 

M 

Rectangle 

58P 

s 

CC 

U.O.(SD) 

M 

Rectangle 

Table  3 


i 

Classifications  of  Problems 

According  to  Smoothness  of  the  Operator  and  Right-Side 

(A=Analytic;  C=Constants;  CC= Computationally  Compli¬ 
cated;  DD=Discontinuous  Derivatives;  E=Entire; 

0= Oscillatory;  P=Peak;  S=Singular) 

Smoothness 
Operator  Right-Side 

Problem  Numbers 

C 

c 

7.  31.  34,  35.  38,  42,  43,  45.  46,  50.  55 

C 

E 

4,  5,  9,  10,  33,  53 

C 

A 

11.  17 

C 

3.  8,41 

C 

S 

3,  14.  47 

C 

0 

6,  11,  53 

C 

P 

10 

E 

E 

12.  21,  37 

E 

A 

2.  6,  18,  40 

E 

S 

54 

A 

C 

16.  27 

A 

E 

1 

A 

A 

26.  30.  32 

A 

S 

15 

DD 

C 

28 

DD 

13.  23,  25 

DD 

S 

22,  25 

S 

C 

29,  51,  56 

S 

s 

19,  24,  36 

o 

A 

6 

|  0 

0 

12 

!  c 

CC 

17.  18 

!  S 

CC 

56 

i  cc 

c 

52 

i  cc 

p 

20 

j  cc 

j 

cc 

30,  39,  44,  48,  49,  54 

Table  4 


Classifications  of  Problems 
According  to  Smoothness  of  the  Solution 


Solution  Smoothness 

Problem  Numbers 

Entire 

1.  4,  5,  9,  10,  12,  19,  20,  21,  22, 

25,  31,  32,  33,  34.  35.  37.  38.  43.  50, 
53 

Analytic 

2,  6.  11.  17,  18,  30,  36,  40.  42 

Singular  Derivatives 

3.  7,  13.  14,  15,  41,  47,  51,  54,  56 

Oscillatory 

6,  11,  12,  33,  35,  38,  42,  50,  53.  56 

Wave  Front 

8,  17,  18.  23.  28.  51 

Discontinuous  Derivatives 

8,  23 

Singular 

54,  55 

Boundary  Layer 

7,  9.  15.  16,  27.  29.  44.  45.  46.  49 

Peak 

10.  20,  24 

Tabled  Solution 

16,  24,  26,  27,  28,  29,  39.  44.  45.  48, 
48,  49,  51,  52 

MM 


APPENDIX  TWO:  MATHEMATICAL  DESCRIPTIONS  ’"I 


PROB  1  Artificial  [7,12,13] 

(e^u  )  +  (e  i  )  -  u/(l  +  x  +  y)  =  f 

xx  y  y 

DOMAIN  unit  square 

BC  u  +  au  =  g 

TRUE  __  xyN  .  .  ...  , 

•  75e  sin  (irxjsin  (ry ) 

Operator:  Self-adjoint,  analytic 
Right  side:  Entire 

Boundary  conditions:  Mixed  except  for  a  =  0. 
Solution:  Entire,  independent  of  a. 

Parameter:  a  introduces  normal  derivative  into 
boundary  conditions. 


PROB  2  Artificial  [12,13] 

u  +  (1  +  y2)u  -  u  -  (1  +  y  )u 
xx  yy  x 

DOMAIN 

unit  square 

BC 

U  -  u  =  q 

N  * 

.135(eX+y  +  (x2  -  x)2log(l  +  y2) ) 

TRUE 

Operator : 

Entire 

1  Right  side:  Analytic 

I  Boundary 

conditions:  Mixed 

!  Solution: 

Analytic 

1  Parameter 

:  None 

PRoeui*  i-  i 


//// 

ih 

S,  1  ! 
\\\ 


.000  .sot 
PR30LEM  2 


PROB  3 


Artificial  1 1 3 i 


u  +  u  =  f 
xx  yy 

DOMAIN  unit  square 
BC  u  =  0 

TRUE  .  a/2  a/2  .  ...a/d-a)  l/d-a),2 

c(x  -  x)  (y  ~  y) ,  c  =  1/  (a  -a  ) 

Operator:  Laplace 

Right  side:  singular  for  a  <_  3 

Boundary  condition:  Dirichlet,  homogeneous 

Parameter:  1  <  a  <_  5  adjusts  singularity  strength 


-f.“  3-  2 


PR8&.EM  3-  3 


5 

n  } 
3*  „  I 
uir-  R  ( 


.’too  TO 

X 


i  .one  to 


PROB  4 


Artificial  [7,12,13] 


u  +  u 
xx  yy 


6xy  e  1 (xy  +  x  +  v  -  3) 


DOMAIN  unit  square 
BC  u  =  0  for  x* 

TRUE  „  x+y 


u  =  0  for  x^O;  u-a(y-y  )u  =g  for  x  =  0 

,  x+y  2  2 .  X 

3e  (x  -  x  )  (y  -  y  ) 


Operator :  Laplace 
Right  side:  Entire 

Boundary  conditions:  Mixed  except  for  a  =  0 
Solution:  Entire,  independent  of  a 
Parameter:  a  introduces  normal  derivative  into 
boundary  conditions 


PROB  5 


Artificial  [13,14] 


4u  +  u  -  au  =  f 
xx  yy 

DOMAIN  unit  square 
BC  u  =  0 

TRUE  2(x2  -  x)  (cos  (2ny)  -  1) 

Operator:  Constant  coefficient,  separable 
Right  side:  Entire 

Boundary  conditions:  Dirichlet,  homogeneous 
Parameter:  a  makes  operator  more  singular  without 
affecting  solution 


PROB  6  Stratospheric  physics  [13,14,16] 

u  +  u  -  (100  +  cos(2trx)  +sin(37ry)]u  =  f 
xx  yy 

DOMAIN  unit  square 
BC  u  =  0 

TRUE  2  4 

-0- 31  (5.4-cos  (4nx) )  sin  (ttx)  (y  -y)  (5.4-cos  (4~y) )  (1/(1+V>  )-.5) 
=  4(x-  .5)2  +  4(y-  .5)2 

Operator:  Entire,  oscillatory,  somewhat  singular 
Right  side:  Analytic 

Boundary  conditions:  Dirichlet,  homogeneous 
Parameter :  None 


PROB  7 


3rtif  i<_  ial 


u  +  u  -  1 
xx  yy 

DOMAIN  unit  square 

BC  u  =  0  _g 

TRUE  Approximate  series  solution  gives  10  '  accuracy 

Operator:  Laplace 

Right  side:  Constant 

Boundary  condition:  Dirichlet,  homogeneous 
Solution:  Has  logarithmic,  singularities  at  corners  in 

second  derivatives;  approximate  solution  is  a 
oolynomial . 


FEATURES' 


iwmiifli . . iniil 


PROB  13 


ntooLEN  13 


( (1  +  (x  -  .4)  u  )  +  u  =  f 

+  x  x  yy 

DOMAIN  unit  square 
BC  u  =  g 

mm[x+*3,  .7+.  5  (x-.4)+ (x-.4)  /(l+x“)  ]  (1+  (y-ll^e  y) 

Operator:  Self-adjoint,  discontinuous  coefficients. 
Right  side:  Line  of  singularities  along  x  =  0.4 
Boundary  conditions:  Dirichlet 
Solution:  Derivative  in  x  is  singular. 

Parameter :  None 


PW8LEH  m-  1 


unit:  square 
u  =  g 

fy“  cos (xv2)  -  l]x2(x  -  l)2 

'•or :  Laplace  plus  nearly  singular  der 
side:  Singularity  in  3-1  v-derin 
gular  for  smal 1 
ary  conditions.-  DirichW 

ion:  B  :  n  _  ridri  i 


elated 


= 0,  aerivati 
ir  sinqulari 
ity  in  solut 


PROB  18 


Artificial 


u  +  (1  +  xv)u  +  cos(x)u  -  e  u  4  3u  = 

XX  VV  XV 


DOMAIN  unit  square 
BC  u  =  g 

TRUE  e-ly2M<.«*>3/U*<3x)3»2l  +  si„(x  .  y  , 

Operator:  Entire 

Right  side:  Large  values  for  x  near  .15 
Boundary  conditions:  Dirichlet 

Solution:  Sharp  wave  front  near  x  =  .15,  entire. 

Solution  if  the  sam.''  as  :r.  the  preceding  problem. 
Parameters:  ■,  id;:'*  *:•»  ••*  renqth  and  shape  of  the 

wave  front. 


. . 


-  PWLB1  1»-  1 

PR0B19  Nonlinear  laminar,  non-Newtonian  flow  [1]  - 1 

(wu  )  +  (wu  )  =  f,  w  =  [T2  +  T2]3  - - - 

X  x  y  y  X  y  8., 

DOMAIN  [-5,1]  x  [.5,1] 

BC  u  =  0  for  x,y  =  1;  u  =0  for  x  =  .5,  u  =0  for  y  =  - 5  «  t 

TRUE  sin(^x)sin(ny)  X  ^  ^N. 

►  '  N. 

Operator:  Self-adjoint,  possibly  singular-  g  , _ 

Right  side:  Possibly  singular.  '"s. 

Boundary  conditions:  Mixed,  homogeneous.  — \ 

Solution:  Entire,  T  similar  to  that  of  non-linear 

problem.  \ 

Parameter:  -.5  <  o  <  1  is  a  physical  parameter.  a  \ 


PROB  20 


From  u  +  u  =  e  [1] 
fcV  xx  yy  x 

u  +  u  -  wu  =  f ,  w=e 
xx  yy 


DOMAIN 

BC 

TRUE 


[0, .51  x  [0,-75] 
u  =  g 


-100(x-.5)‘ 


10<p  (x)v  (y)  +  a  where  <p  (x)  =  e  ( 

Operator:  Helmholtz  type,  approximates  nonlinear 
operator. 

Right  side:  Sharp,  large  values  near  x  =  y  =  -5. 
Boundary  conditions:  Dirichlet,  homogeneous. 
Solution:  T  lias  a  peak  at  x  =  V  =  .5. 

Parameter:  a  adjusts  singularity  of  operator. 


(x  -x) 


PROB  21 


HWXJE*  21 


Artificial 


Au  +  Bu  +Cu  =f,A=C=l+ 
xx  xv  yy 
DOMAIN  unit  square 

BC  u  =  g 


TRUE 


u  =  g 
x+y 


Operator:  Entire,  has  mixed  derivative  term 

Right  side:  Entire 

Boundary  conditions :  Dirichlet 

Solution:  T  is  entire 

Parameter :  None 


PROB  22  Elastic-plastic 


torsion  [15] 


w(u  +u  )+wu  +  wu  =  f,  w  defined  below 
xx  yy  xx  y  y 

DOMAIN  unit  square 
BC  u  =  g 

TRUE  [17.06  +  3.62  (x2  +  y2))  (x2  -  1)  (v2  -  1) 

Operator:  Expanded  form  of  self-adjust  problem-  discon¬ 
tinuous  coefficients,  w  =  1/7996  if  A  .0025 

w  =  1/(236  +  19.4/A)  if  A  >  .0025  where  A  =  A2  *  rj.  r; 
Right  side:  Singular  " 

Boundary  conditions:  Dirichlet 
Solution:  7  is  a  quartic  polynomial 
Parameter:  None 


PROB 


Friction  in  a  brak ■■  sr.oe  { 3 ] 


PROBLEM  24-  1 

^  (R=0.2S,  8=100 ■ .  C=-0,10) 


u  +  u  + — -  u  =  vh  /h3,  h  =  sin(anxy) 
xx  yy  h  x  x 

DOMAIN  [-1,1]  x  [-1/H 
BC  6u  +  u  +  u  =  0 


TRUE 


Unknown 


Operator:  Laplacian  plus  u^_  term  which  is  possibly 
singular . 

Right  side:  Analytic,  possibly  nearly  singular. 
Boundary  conditions:  Mixed,  homogeneous. 

Solution:  Approximate  solutions  given  for  8  cases. 
Parameters:  a, 6  and  y  are  physical  parameters, 


1.  a  =  .25  6  =  100  Y  =  -.1 

2.  a  =  .25  6  =  1000  Y  =  *1 

3.  a  =  .5  8=1  Y  =  --01 

4.  a  =  .5  6  =  10  Y  =  --1 

5.  a  =  1  6=1  Y  =  "-I 

6.  a  =  1  6  =  10  Y  =  “-1 

7.  a  =  1  6  =  100  Y  =  -1- 

8.  a  =  1 


6  =  1000  Y  =  -1' 


PROBLEM  24-  2 

1  (fl=0.2S.  B=1000..  C=-0.10l 


PROBLEM  24-  3 
|  !  H=0 .50 »  8=1.,  C=-0.01) 


PROBLEM  24-  4 
8  (RrO.SO.  6=10..  C=-0.t0) 


PROBLEM  24-  6 
I  lffcl.00.  8=1..  C=-0.10) 


PROB  25 


Artificial 


a  a  a-1  a-1  ,  ,  .a  _ 

-xu  -  y  u  -ax  u-ay  u+  (xy)  u  =  f 
xx  yy  x  y 

DOMAIN  unit  square 
BC  u  =  0 

TRUE  3ex+y(x  -  x2)(y  -  y2) 

Operator:  Variable  siroothness,  exoanded  self-adjoint 
Right  side:  Variable  smoothness 
Boundary  conditions:  Dirichlet,  homogeneous 
Solution:  Entire,  does  not  depend  cn  parameter  a 
parameter:  a  affects  smoothness  of  operator  and  right 
side  without  affectina  solution. 


PROBLEM  25-  I 
|  lfel.5) 


PROB  26  Viscous  flow  [31 


2  3  2  ^ 

uxx+u  +Aux=-60ax/B  where  B=  (a+x  )  ,  A=6x(l+x  )  /B 

DOMAIN  10, a]  *  [0,a] 

BC  u  =  0 

TRUE  unknown 

Operator:  Laplacian  plus  u  term.  For  a  =  1  it  is 
expansion  of  a  self-adjoint  operator. 

Right  side:  Analytic 

Boundary  conditions:  Dirichlet,  homogeneous 

Solution:  Approximate  solutions  found  for  a  =  1,5 
and  10. 

Parameter:  a  is  a  physical  parameter  adjusting  the 
domain  and  entering  the  coefficients. 


PROBLEM  26-  1 
(HzJ) 


PROBLEM  26-  2 
(R:S) 


PROBLEM  26-  3 
8  IfizlO) 


PROB  27,  Distribution  of  diffused  par'  i'-lc-s  [3[ 

u  +  —  u  +  -or  u  +  -or  (cot  y)  *u  =  -100 
XX  x  x  2  yy  _2  y 

DOMAIN  [.1,1]  x  (.1,1] 

BC  u  =  0 

TRUE  unknown 

Operator:  Nearly  singular,  analytic 
Right  side:  Constant 

Boundary  conditions:  Dirichlet,  homogeneous 
Solution-  Annroximate  '■!  i o- **'-> . 

Paran-'f-1-  t-no 


PROB  28  Artificial 


(wu  )  +  (wu  )  =1  where  w  =  a  if  0  £  x,y  1  -5 

x'x  H  =1  otherwise 

DOMAIN  [-1,11  x  t-1/1] 

BC  u  =  0 

TRUE  unknown 

operator:  Self-adjoint,  discontinuous  coefficients. 

Right  side:  Constant 

Boundary  conditions:  Dirichlet,  homogeneous 
Solution:  Approximate  solutions  given  for  a  =  1,  10, 
100.  Strong  wave  fronts  for  a  »  1. 

Parameter:  a  adjusts  size  of  discontinuity  in  operator 

coefficients  which  introduces  large,  sharp  jumps  in 
solution . 

PRO0LEH  28-  2  PROBLEM  28-  3 


PROB  29 


uao 


Many  physical  interpretations  [10] 
.  a 


PROBLEM  29-  2 


PROB  4  30 


Artificial 

1 


[2+  (y-1) e  y  lu  +tl+ - ■ 

XX  l+(2x)P  Yy 
DOMAIN  unit  square 
BC  u  =  g 

TRUE  „  ^  .2  * 


]  u  +y [x (x-1)  +  (y- . 3) (6-.7)]u 

yy  =f 


x  j-  y 


+  (y-1) (l+x)e 


+  Y (x+y) cos (xy) 


1+ (2x) 

Operator:  Coefficients  may  be  widely  varying,  singular. 
Right  side :  Complicated  behavior 
Boundary  conditions :  Dirichlet 

Solution:  Complicated  behavior,  with  wave  fronts,  etc. 
Parameters:  a,  6,  y  adjust  the  contribution  of  3 
independent  complexities  of  the  problem. 


PROBLEM  30-  1 
|  (ffcl.O.  8=6,0.  Cs-1.0) 


PROBLEM  30-  2 
8  (R=I.O.  B=3.0.  C=0 .51 


PROBLEM  30-  3 
8  IR=23.0.  6=2.0.  C 


PROBLEM  30-  4 


EHiawaarar 


PROBLEM  30-  S 
8  IfclOO.O.  6=5.0.  C=2.0J. 


PROBLEM  30-  6 
8  (fclO.O.  B=M.O.  Csl.O) 


PROBLEM  30-7  ,, 

(IfcM.O.  B=S,0.  Cr-O.S) 


PROBLEM  30-  B 
|  (R=3.Q.  B=6.0»  C=2.0) 


PROBLEM  30-  9 
8  lfl=0.S.  B=3.0.  C=)0.0) 


.'.2-1  3 


PROB  31 


Temperature  distribution  [5] 


u  +  u  =  -1 
xx  yy 

DOMAIN  [-1,1]  *  [-1,1]  S 

BC  u  +  u^  =  g 

TRUE_  (x2+y2)/4  +  .Q21564  +  . 01440 (x4-6x2y2+y4)  ® 

+  .0000493  (x8-28x6y2+70xV-28x2y6+y8)  *g 

-.00000064  (x12-66x1°y2+495xV-924x6y6+495x4y8-66x2y1  ‘‘ 

Operator :  Laplace  . 

Right  side:  Constant 
Boundary  conditions:  Mixed 

Solution:  Harmonic  poly,  expansion  for  homo.  BC.  5 

Parameter :  None  51 


PROB  32 


32  Stress  in  helical  spring  [5] 

u  +  u  +  — - —  u  =  f 
xx  yy  5  -  y  y 


-  PROBLEM  32 

*1  MRX:  .127H  MIN:  .0000 


DOMAIN  [-.5, .5]  x  [-1,1] 

BC  u  =  0 

TRUE  (i-y2)(i-4x2)  (5-y3)  (.0004838y  +  .0010185) 

Operator:  Analytic 
Right  side:  Analytic 

Boundary  conditions:  Dirichlet,  homogeneous 
Solution:  Polynomial  obtained  by  Ritz  method  for  a 
physical  problem. 

Parameter:  The  5  in  the  operator  is  a  value  of  a 
physical  parameter. 

PROB  33  Torsion  on  a  shaft  [5] 

u  +  u  =  f 
xx  yy 

DOMAIN  [0,1]  x  [-1,1] 

BC  u  =  g _ 

TRUE  P  =  14  +  ^133,  q  =  14  -  /L33,  r  =  (7-q)/ (r*533) , 

t(y)  =  1-y2,  C  (x)  =  e  pX  -  e  qX,  B(x)  =  (7-p)r/16C  (x) , 

A  (x)  =  rC(x)  +  e  qX,  TRUE  =  t(y)[A(x)  +  t(y)B(x)] 
Operator:  Laplace 
Right  side:  Entire 
Boundary  conditions:  Dirichlet 
Solution:  Entire 


PROB  34  From  infinite  region  problem  [5] 

u  +  u  =  -1 
xx  yy 

DOMAIN  1-1.1]  x  1-1,1] 

BC  u  =  g 

TRUE  .295776  -  (x2+y2 )/4  -  14476  (x4-6x2y2+y4)/319424 

+  429 (x8-28x6y2+70x4y4-28x2y6+y8)/319424 
Operators :  Laplace 
Right  side:  Constant 
Boundary  conditions:  Dirichlet 
Solution:  Harmonic  polynomial  expansion  for 

homoqenr-o  i~.  boundary  conditions. 

Parar*  *  r  ■ 


MRX:  T.292S  MlN=  -1.1169 


X  [-1 
y  =  ±1 
80  Ip  + 
2  4 


PROB  36 


Adapted  from  Problem  27 

2.1  .  cot  y  .. 

—  u  + - —  u  + -  u 


(1  +  8)u  + -  u  + - -  u  + - -  u  =  f 

1  xx  x+a  x  a)2  yy  x+a  y 
DOMAIN  unit  square 
BC  u  =  g 

TRUE  (1  -  B)eX+y  +  Blogp(x  +  a) 

Operator:  Possibly  singular  coefficients  for  a  =  0. 
Right  side:  Analytic  except  for  a  =  0;  then  singular. 
Boundary  conditions:  Dirichlet 
Solution:  Logarithmic  singularity  for  a  =  0. 

Parameters:  rt  adjusts  distance  of  singularity  from 

domain,  3  adjusts  relative  size  of  exponential  and 
logarithmic  terms  in  solution. 
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PROB  37 


From  nonlinear  minimal  surface  [1] 


Au  +  Bu  +  Cu  =f,  A=  (1+T  )‘ 
xx  xy  yy  y 

DOMAIN  unit  square 

BC  u  =  g 

TRUE  T(x,y)  =  (x  -  3y)2e*  y 

Operator:  Includes  cross  derivative. 

Right  side :  Uitire 
Boundary  conditions:  Dirichlet 
Solution :  Entire 

Parameter:  None 


B=-2TxTy,  C  =  (1+Tx)‘ 


PROB  38 


Electrostatics  [11] 


u  +  u  =0 
xx  yy 


DOMAIN  t-n/2*ir/21  x  [0,1] 

BC  u  =  g  for  x  =  ±tt/2,  y  =  1;  u^  =  g  for  y  =  0 

TRUE 

e  cos [(2a+l)x]sinh[(2a+l)y]/ (2a+l) 


Operator:  Laplace,  homogeneous 

Right  side:  Zero 

Boundary  conditions :  Mixed 

Solution:  Entire,  may  be  oscillatory. 

parameter:  a  adjusts  the  oscillations. 


,  pR0«.EN  36-  3 

|  FEWuRES:  000.26  000.00  000.60  015.60 
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PROB  41  -  [2oi 


u  +  u 
xx  yy 


f  >u  =  I 


DOMAIN  [0,n]  *  [0,ir] 

BC  u  =  0 

TRUE  approximate  solution  accuracy  depends  on  8 
x(n-x>  4  £  sin  [ (2k-l)x]cosh  [  (2k-l)  (y-ii/2)  1 

k=l  (2k-l)  '"cosh  [  (2k-l)iT/2] 

Ope  rat or :  He lmhol tz 

Right  side:  Series  for  function  with  singularities. 
Boundary  conditions:  Dirichlet,  homogeneous. 

Solution:  Infinite  series  converging  like 

1/k^.  The  solution  has  derivative  sinqularities. 
Parameters:  a  adjust  u  term,  possibly  makes  operator 
nearly  singular.  8  is  number  of  terms  in  series. 


PROB  43  Artificial  [20] 


u  +  u  +  u  =  0 
xx  yy  x 

DOMAIN  [0,-]  x  [o,u] 

BC  u  =  0  for  x  =  0,  y  =  0 , n ;  u-u  = ‘2  sin  (y+n/4 )  for  x 


TRUE 


e(n-x)/2  sinh  (,'5x/2)siny 
sinh \&-z/2 ) 


Operator:  Constant  coefficient,  homogeneous. 

Right  side:  Zero 

Boundary  conditions:  Mixed 

Solution:  Entire 

Parameters :  None 


PROB  44 


From  nonlinear  problem  i2u] 
u  +  wu  =  w 


DOMAIN 

BC 

TRUE 


unit  square 
u  =  0 
unknown 


Operator:  Helmholtz  type 

Right  side:  Complicated 

Boundary  conditions:  Dirichlet,  homogeneous 
Solution:  Approximate  solution  given  for  r  =  r(x,y) 
tabulated  from  a  solution  to  the  nonlinear  problem: 
should  be  u, 

,  ,  2  n  .3-1  [ySr/d+yr)] 

w  (x,y )  =  -a  (1-r)  e 

Parameters :  a,S,  V  and  5  are  physical  parameters. 


w (x,y )  =  -a  (1-r)  e 
Parameters:  a,S,  Y  anc 
Four  cases  are  given: 

1.  a  =  1-425  3=1 

2.  a  =  10  3=1 

3.  a  =  1.425  3=2 

4.  a  =  1.425  3=2 


’S33.*"  — -  ? 
•irlO.Ci  3=1  ■ 


are  physical  parameters. 

=  .5  5  =  2 
=  .5  5=2 
=  ,~4  5  =  25 
=  .5  5  =  2 


I  0=2.  Cz.&t.  Ct2S; 


-- 

I  :fi=;.«i25.  i 
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PROB  47  Artificiai 

u  +  u  =  f 

s  , 

xx  yy 

DOMAIN 

unit  square 

BC 

u  =  g 

i. 

TRUE 

/  ^/2 
(xy) 

Operator : 

Laplace 

§j 

Right  side 

:  Variable  singularities 

I  Boundary  conditions :  Dirichlet 

£ ' 

Solution: 

Singularity  of  variable  strength  - 

Parameter: 

a  adjusts  singularity  strength. 

B;  .  . 

.k  -23C 

PR3E.EM  2 

PROBLEM  47-  3 

PRS-E**  46 

PROB  48  Nonlinear  diffusion  in  catalysts  [2j 

,  (Y-l)  (as  (1-r) )/ (1+3  (1-r))  _  n 

u  +  u  -1.425r  e  u  -  0 

xx  yy 

DOMAIN  unit  square 
BC  u  =  1 

TRUE  Unknown 

Operator:  Helmholtz  type,  homogeneous 

Right  side:  Zero 

Boundary  conditions:  Di rich let 

Solution:  Approximate  solution  given  for  r  =  r(x,y) 
tabulated  from  a  nonlinear  PDE  solver;  r  should  be 
l  Parameters:  (a,3,v)  are  physical  parameters.  5  cases 
given:  (1,-04, 2),  (1,-04,25),  (2, .04, 2),  (2,-04,25) 

and  (2, .5, 2). 


>9-  3 

i  i«;2.  8: .04.  C=2I 


PROBLEM  H6-  R 
I  ife2S.  fc.S.  fc2J 


"“SSlE!*  -S-  2 
!f=25-  Sr.C-. 


(  ■  /  *  « 
6,  ' 


u-e,  v 
6.  ■ 


-9*  S 
B  .S3.  : 


occ  zx  ,*oc 


PROB  51 


b  "1=0.3; 

e - r— R-r-f- 


U  “  U  t  ■  „  -  -10 

XX  XX  2  yy 

DOMAIN  unit  square  x 

BC  u  =  0  for  x  =  1;  u  =  0  for  x,y  =  0;  Au  +Bu  =  0  for  y  = 1 
TRUE  Unknown  T 

Operator:  Singular  coefficients 
Right  side:  Constant 
Boundary  conditions:  Mixed 

Mxl  -  (?  *  >  ”  B(x)  .  *  >  » 

(1  X  £  a  [0  X  <_  a 

Solution:  Has  singularity,  unusual  behavior. 

Parameters:  a  adjusts  position  of  change  in  boundary 

condition  for  y  =  1. 

PROBLEM  51-  2  PROBLEM  51-  3 


PROBLEM  51-  M 


\  \  \  ! 


PROB  52 


Nonlinear  reaction  [2] 


r(u  +u  )  +•  r,  u  +  r  u  -  au  =  0 
xx  yy  lx  2  y 


DOMAIN  unit  square 

BC  u  +  u  =1 

TRUE  _ _ 


Unknown 

Operator:  Expanded  from  self-adjoint  PDE,  homogeneous. 

Right  side:  Zero 

Boundary  conditions :  Mixed 

Solution:  Approximate  solution  for  r(x,y)  tabulated 

[  from  nonlinear  PDE  solver;  r  should  be  11/(1 +  lOu) ; 

r,  ,r„  are  finite  differences  for  r  ,r  . 

12  y  y 

Parameter:  a  =  d1.  a!  •.aramc  *  <  r . 


PROBLEM  52-  1 


PROBLEM  52-  3 


PROB  53 


Artificial 


xx  yy 

DOMAIN  unit  square 

BC  u  =  g 

TRUE  cos  (By)  sin (6 (x-y) ) 

Operator:  Helmholtz 

Right  side:  Entire 
Boundary  conditions:  Dirichlet 
Solution:  Entire,  oscillatory 

Parameters:  a  can  make  operator  nearly  singular-  p 
adjusts  the  oscillations  of  the  solution. 


PROBLEM  63-  2 
(fiUO.O.  BrPI) 


PROBLEM  S3-  3 


PROB  54  Artificial 


(1+X2)u  +(1+A2)U  +2xu  +16yAu  - (1+ (8y-x-4)  )u  =  f 

xx  yy  x  y  A(y)  „  4y2  +  0 

DOMAIN  unit  square 

BC  u  =  g 

TRUE  a  =  max[0,  (3-x/A(y) )  ],  C  =  max [0,x-A (y)  ] 

D  =  0  if  C  ^  .02 ,  D  =  e  B^C  if  C  ^  .02 
u(x,y)  =  2 .25x (x-A (y) ) 2 (1-D) / (4A (y) 3 ) +1/ (1+ (8y-x-4 )  ) 
Operator:  Expanded  form  of  self-adjoint  operator. 

Analytic. 

Right  side:  Complicated  with  possible  wild  behavior. 
Boundary  conditions:  Dirichlet 

Solution:  Wildly  behaving  for  a  possible,  has 
singularities  for  x  -  4y2  =  a  or  4y  =  -a. 


ss- 
(A 
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PROB  55  Conducting  fluid  in  magnetic  fir  id  :  ■ 

u  +  u  =0 
xx  yy 

DOMAIN  [0,6]  x  [0,1] 

BC  u  =0  for  x=0;  u=f  for  x=6;  Au+Bu  =g  for  v=0, 

TRUE*  Unknown  ^  Cu+Du^=h  fir  y=l 

Operator:  Laplace,  homogeneous 
Right  side:  Zero 

Boundary  conditions:  Mixed  and  complicated: 
f  (y )  =  ^'/^sin(ny/2)/a^ 


fl  x  <  a 
\x  x  a 


B  (x)  = 


C  (x)  = 


x  <  a  or  k  =  1 
x  >  a  and  k  =  2 


D  (x)  = 


0  x  <  a  or  k  -  i 
1  x  >  a  and  k  =  2 


=  / 

I  2  (a-x 


^/a3  x>a 


H(x)  =  < 


jO  x<a  or  k=l 

\  3(ct-x)  3 

ve  /a  \>a  and  k=2 


Solution:  das  singularities  at  boundaries,  widely  vary 
ing  behavior. 

Parameters:  a.  8  are  physical  parameters,  c  selects 
different  ohysical  models.  Four  cases  are  Given. 


‘■RCBiEt'  55-  3 
Iftel..  B:3  .  C=Si 

FEST5RES:  308.15  CCO.OC  COS  2C  C2".'< 


PROB  56 


Art:  la'. 


1  1 

u  +  —  u  +  —  u  =0 
xx  x  x  2  yy 
x 


xx  x  x  -i  yy 
DOMAIN  [0,1]  x  [o,2jt]  * 

BC  u  =  g 

TRUE  a  -z  xjiny  -z^xcosy  7  2n 

l  w  [e  cos  (z.xcosy) +e  K  cos (zkxsiny / ,  ■/.  a^x  cos(2ny) 

k=0  k  -1 

Gauss  weights  w  and  points  z^  depend  on  a,  a^  =  (1,-5,  ✓  i/6, -1/10, -1/15, 1/30,-1/50) 

Operator:  Singular  coefficients,  homogeneous* 

Right  side:  Constant 
Boundary  conditions:  Mixed 

Solution:  Series  expansion  approximates  electrosta- tcs  solution. 

Parameters:  a  =  order  of  Gauss  quadrature  for  integral,  $=  no.  terms  in  expansion. 
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